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Abstract 

We show that universality in chaotic elements can be 
lifted to that in complex systems. We construct a glob- 
ally coupled Flow lattice (GCFL), an analog of GCML 
of Maps. We find that Duffing GCFL shows the same 
behavior with GCML; population ratio between syn- 
chronizing clusters acts as a bifurcation parameter. 
Lorenz GCFL exhibits interesting two quasi-clusters 
in an opposite phase motion. Each of them looks like 
Will o' the wisp; they dance around in opposite phase. 
Keyword : synchronization, universality, globally cou- 
pled maps and flows, Lorenz model, Duffing oscillator 

1 Introduction 

In studying chaotic systems we may consider two facts 
as guiding principles. One is 'Universality in Chaos' pQ . 
Especially we here consider the universality between a 
map (e.g. a logistic map, a circle map • • • ) and a 
flow, that is a system described by an ordinary differ- 
ential equation (a Duffing oscillator, Lorenz flow • • • ). 
The other is a self-organization of a certain attractor 
by synchronization 2 in a complex system under con- 
flict between randomness and coherence. One impres- 
sive model embodying this phenomenon is Globally 
Coupled Map Lattice (GCML) [3], g|. In this paper 
we study Globally Coupled Flow Lattice (GCFL) [5] 
and show that the universality at the level of the con- 
stituents may be extended to the level of the whole 
system, that is, our GCFL has much in common with 
GCML. 



2 Universality in Chaos 

Let us first briefly recapitulate the universality in 
chaos at the level of elements pQ. Perhaps the sim- 
plest chaotic system is an iterated logistic map £ n +i = 
fa(%n) with a quadratic map f a (x) = 1 — ax 2 and a 
nonlinear parameter a in the range < a < 2. An- 
other simple system is an iterated circle map where the 
map is replaced by gb(x) = 6sin(7nr). Both systems 
share the same Feigenbaum ratio 4.699 • • • and subject 
to the same universality class. This is a piece of the 
universality in chaos, which was beautifully explained 



by Feigenbaum and by Cvitanovic using renormaliza- 
tion group [1], [6]. 

Next let us look at a system described by an or- 
dinary differential equation (a flow for brevity). The 
Lorenz model, for instance, is described by a three di- 
mensional ODE 



dx 
~dt 



a(y - x), 



dy 
dt 



-xz ■ 



dz 
~dt 



xy — bz. (1) 



We here consider the range of r from 230 to 200 
(cr = 10, b = |), where the attractor bifurcates from 
a limit cycle to chaos. As for the 2-dimensional flow, 
we consider a Duffing oscillator, 



dx dy 



■ X 



Acos(t). 



(2) 
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chaos. Above three models exhibit the 
furcation structure as well as periodic 
The reason why is essential 
The attractor of Lorenz model 
reduces to two dimensional due 
to dissipation. Then, the 
Poincare section of Lorenz flow 
(that of Duffing oscillator) is 
nothing but one-dimensional 
iterated map, see Fig. [TJ That 
is, if the one-dimensional map 
bifurcates, the corresponding 
flow must also bifurcates. And Figure 1: 

then, at the bifurcation limit, one-dimensional maps 
(in the same class) are all governed by a universal 
function. Therefore, we can amazingly understand the 
route from convection to turbulence from a simple lo- 
gistic map. 



3 GCML 

The GCML has N maps as elements and evolves under 
all to all interaction via their mean field. Here we 
consider its simplest form, where the map is a logistic 
map with the same a and the coupling is also common 
(homogeneous GCML) [3] @]. It evolves in an iteration 



of two steps. The first is a parallel mapping 

xf d = f a (xi(n)) = 1 - ax 2 t (n), i=l,---,N (3) 

where each map evolves with high nonlinearity and 
randomness is introduced in the system. The mean 
field is then calculated as h(n) = Y%=i xf d /N. The 
second step is interaction via the mean field 

Xi (n + 1) = (1 - e)xf d + eh{n), i = 1, • • • , N (4) 

with a coupling constant e. In this step all the maps 
are pulled to the mean field and coherence is in- 
troduced in the system. Under the conflict of ran- 




Figure 2: 

domness and coherence GCML shows various inter- 
esting phases. The phase diagram, explored first by 
Kaneko[3j is shown in Fig. [3j We should note that 
the region with a very small e is called the turbulent 
regime, but actually there occur drastic global periodic 
motions of maps, if e takes certain values for a given 
a. The turbulent regime is controlled by the periodic 
windows and their foliations [3]. 

4 Duffing GCFL 

We transport the basic structure of GCML, the iter- 
ation of independent evolution and subsequent inter- 
action, to GCFL. We choose as elements N Duffing 
oscillators that are discretized in time. The first step 
is 

x%f=x {i) {t)+y {i) {t)At 

y™) d = V(i)(t) + (-ky { i)(t) - (*(;)(*)) 3 + 7.5cos(t)) At 

(5) 

and the second is 

N 

X(i) {t + At) = {l-e)x%f + £ fY, x tf 

1=1 (6) 

N v ' 

y{i) (t + At) = {l-e)y™? + £ fY { y™f. 

In order to investigate the phase structure of Duffing 
GCFL, we need some rough estimate of the starting 



point in the parameter space, since otherwise Ed, for 
instance, can be any value from to 1. In GCML the 
two clustered phase in opposite phase is the most re- 
markable state formed by synchronization and realized 
for e « 0.19 — 0.27 at a « 1.68. Let us start estimat- 
ing corresponding values of Duffing parameters k and 
Ed- Whether corresponding cluster state is realized in 
Duffing GCFL or not is of course highly nontrivial. 

A natural way to estimate k is to use the univer- 
sality at the element level and match the bifurcation 
tree of the Duffing oscillator with that of the logistic 
map by a scale transformation^- We find that Duffing 
k « 0.24 corresponds to logistic a w 1.68. 

As for ed, we should notice that the averaging in- 
teraction is introduced at every At in GCFL. Hence 
,if the period is T, the estimate may be given by 

(l_ £D )23t =l-e map . (7) 

The inclusion of factor 2 in the denominator needs 
short explanation. The two clusters in GCML are os- 
cillating in opposite phase and at every two steps each 
comes back to previous value. 

(+- + - + -•••) — + 
(- + - + - + •••) — - 

Therefore two steps (one step) in GCML correspond 
to one period (half period) in Duffing GCFL. With 
T = 2tt and At = 10 -3 , we obtain 

6.6 x 10" 5 < £ flow < 9.7 x 10" 5 . (8) 

We should add that we have approxi- 
mated the sequence of two steps, Eq.© and 
Eq.(j6]), namely (evolution • interaction) T ^ At by 
(evolution) T / At . (interaction) 7 ^ A£ in estimating ed- 

4.1 Phase Structures 

In Fig. [3] we show the phase structure of Duffing GCFL 
with respect to e at A = 7.5 and k = 0.24. Let us 
explain the observed states in order. 

(1) Turbulent phase. At very small £, each flow 
evolves almost independently in a chaotic orbit. 

(2) Many clusters in chaos. Flows divide into many 
clusters and each cluster evolves chaotically. 

(3) Two chaotic clusters. The same as above but the 
number of clusters is only two. 

(4) Two periodic clusters. Flows divide into two clus- 
ters and each cluster evolves periodically in ap- 
proximately opposite phase to the other. 

1 Universality in bifurcation is strictly a valid concept at the 
limit of sequential bifurcation, but in practice a simple scaling 
makes the two trees overlap each other very well. 
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Figure 3: Upper: Phase diagram Duffing GCFL with 
respect to e at A = 7.5, fc = 0.24. Lower: attractors of 
he model (N = 20) in three typical region; turbulent, 
two clustered state, coherent chaos with increasing e. 



(5) Two periodic clusters or coherent chaos depend- 
ing on the initial values of flows 0. 

(6) Coherent chaos. All flows at strong coupling 
bunch together and move in a chaotic orbit. 

To our knowledge, this is the first time observation 
of the fourth phase (dark shaded region in Fig. |3J)— 
periodic clusters of flows in opposite phase motion just 
in accord with GCML [7j. We further study this phase 
below. But not only this but also the global phase 
structure is in agreement with GCML. This is a case 
that universality between elements (a logistic map ~ a 
Duffing oscillator) also holds between systems (logistic 
GCML - Duffing GCFL). Note that the range for the 
two clustered phase is 



9.2 x 10~ 5 < e D < 12.8 x 10~ 5 



(9) 



as seen in Fig. [3l This overlaps with (jHj) but is shifted 
to the larger. Considering that our estimate is rather 
rough, this agreement is remarkable. 

4.2 Two Clustered Phase 

In GCML this phase is a consequence of reduction of 
nonlinearity due to the averaging interaction [H [9l [10] 
and the fluctuation of the mean field is minimized for 
the stability of the system. Further interesting prop- 
erty of the two clustered attractor in GCML is that 
it may be controlled by changing the population ratio 
of two clusters [3l [8]. Let us call N+ (iV_) the number 
of maps beyond the mean field at some iteration step 
(say n = 1000) after random start. Then the attractor 
of GCML is determined if we give not only a and e but 
also = N+/N. 



2 This is because Duffing oscillator has two possible attrac- 
tors mutually symmetric with respect to the origin of the xy 
plane. In the coherent chaos case, all flows are attracted to 
one of the attractors, while in the 2 clusters case, flows divide 
themselves into two groups, one for each attractor. 




Figure 4: Bifurcation of final attractors in two clus- 
tered phase. Upper: Logistic GCML with a = 
1.68, e = 0.23 and N = 10 4 . (a) Evolution of two 
clusters (0 = 0.41). (b) bifurcation of cluster attrac- 
tor with respect to 0. Maps are plotted for even n 
(n = 1000, 1002, • • •)• Lower: (c) Duffing GCFL with 
A = 7.5, k = 0.24, e = 1.1 x 10~ 4 and N = 10 4 . 
(0 = 0.22) (d) The same with (b) but for GCFL. 
Poincare shot of flows are taken at every At = 2tt. 



In Fig. [H(a) we show how the GCML two clusters 
evolve when they are formed with fixed from random 
start and in (b) how their orbits change with by 
plotting all maps Xi(n) at every even n in order to 
separate the two clusters. Correspondingly we show 
in (c) the evolution of GCFL clusters with fixed and 
in (d) the bifurcation structure of the cluster orbits by 
plotting the Poincare shot of flows at every 2tt. We see 
clearly that both GCML and GCFL cluster attractors 
share the same bifurcation structure with the variation 
of 0. In more detail, GCML at even iteration steps 
and Duffing GCFL at every At = 2tt are in one to one 
correspondence. 

5 Lorenz GCFL 

The Lorenz GCFL is constructed just in the same way 
with Duffing GCML (the iteration of two step process, 
discretized evolution in At and subsequent interaction 
via the mean field). The model parameters are now 
r and Ed with a = 10 and 6 = 8/3 fixed. We are 
interested in the self-formation of cluster structures so 
we use large TV, typically N = 10 4 . 

By matching the Lorenz tree with the logistic tree 
we estimate that r = 208 for Lorenz flow corre- 
sponds to logistic a = 1.68. At r = 208, the one 
turn of Lorenz chaotic attractor is roughly T = 0.50. 
With At = 10 -3 , the correspondence gives £l ~ 
8.2 x 10" 4 - 12.0 x 10" 4 . It turns out that there 
occur interesting quasi-periodic cluster attractors for 



£l both at about 50 percent and about twice of this 
range, while at £l within this range there occurs tran- 
sitive behavior between the two. For further small £l, 
flows evolve almost randomly, while for further large 
£l, all flows bunch together and show coherent chaos 
with the high nonlinearity of the element at r = 208. 
This is just the same way with the GCML and Duffing 
GCFL. We stress that only with the help of our esti- 
mate from universality, we can access these interesting 
three regimes. 

Let us investigate the strong coupling regime in de- 
tail!. Here we find interesting behavior of quasi-spatio 
clusters; each of them, if we may say, resembles will- 
o'-the-wisp very much. In Fig. [5] we show the final 
stages of their formation process (projecting into the 
xy plane) from a random start at the central sl of this 
regime (« 24 x 10" 4 ). At around (a) (t « 3) most of 



(b) 
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Figure 5: Formation of clusters in Lorenz GCFL with 
e = 24 x 10- 4 , r = 208, N = 10000, At = 10~ 3 . 
(a) t = 3, (b) t = 30. 



flows divide themselves into two large spatial clusters 
C and C and in addition two small clusters s and s 
with high population density are formed. In the xy 
projection C and C (s and s) locate oppositely each 
other with respect to the origin reflecting the symme- 
try of Lorenz flow. They circulate around periodically 
in the first and third quadrant respectively (always op- 
posite in phase). Hence, the fluctuation of the mean 
field in the xy direction is suppressed. The clusters 
C and C actually exchange some part of their mem- 
bers at every time when they come close each other in 
circulation. Always the spatial clusters C and C exist 
but members are gradually mixed up. In this sense C 
and C may be called as quasi-clusters. (On the other 
hand cluster s and s do not mix.) Around t = 25 (after 
about 100 circulations in the respective quadrant) the 
smaller cluster s (and s) disappear. The quasi-clusters 
C and C remain circulating opposite in phase. 



6 Conclusion 

We have constructed GCFL, matching the nonlinear- 
ity parameter and the coupling with those of GCML. 
The matching is only a necessary condition for the 
GCFL to inherit the properties of GCML. We sum- 
marize below to what extent the intriguing features of 
GCML are realized in the GCFL. 

Duffing GCFL preserves all of features of GCML, 
except some shift of the coupling sd to the larger side 
than the prediction. The final attractor is remark- 
ably controlled by the population ratio 6 just in the 
same way with GCML. On the other hand the flows of 
Lorenz GCFL do move on a two dimensional surface 
but they do not form a tightly bound clusters. They 
form quasi-clusters (with mixing) in opposite phase 
motion at the strong coupling regime. 
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